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ABSTRACT 

We consider the dynamical evolution of an interstellar cloud that is initially 
in thermal equilibrium in the warm phase and is then subjected to a sudden 
increase in the pressure of its surroundings. We find that if the initial plasma (3 
of the cloud is of order unity, then there is a considerable period during which 
the material in the cloud both has a small (3 and is in the thermally unstable 
temperature range. These conditions are not only consistent with observations of 
star forming regions, but are also ideally suited to the production of density inho- 
mogeneities by magnetohydrodynamic waves. The end result should be a cloud 
whose size and average density are typical of Giant Molecular Clouds (GMCs) 
and which contains denser regions whose densities are in the range inferred for 
the translucent clumps in GMCs. 

Subject headings: galaxies:star fomation - MHD - fast shocks 



1. Introduction 

Observations of magnetic fields in star 
forming regions (e.g. Crutcher 1991) tell 
us that the plasma (3 is of order 0.04 and 
that the velocity dispersion is of the or- 
der of the Alfven speed. For this rea- 
son, many of the simulations of the early 
stages of star formation (e.g. Fiedler & 
Mouschovias 1993; Ballesteros-Paredes & 
Mac-Low 2002; Falle & Hartquist 2002; 
Padoan & Nordlund 2002; Gammie et al. 
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2003) have assumed that the initial value 
of f3 is low. These papers show that dense 
structures can be formed under these con- 
ditions, but they do not consider how f3 
came to be so low in the first place. 

On scales large compared to the separa- 
tion between clouds, the thermal and mag- 
netic pressures in the interstellar medium 
are comparable. Locally (3 may be about 
0.5, while in the midplane at a distance of 
5 kpc from the Galactic Centre, it may be 
about 1.5 (Ferriere 1998). The production 
of a star forming region must therefore also 
reduce (3 by a considerable factor. Here we 
examine a mechanism for this which relies 
upon pressure driven compression and ra- 
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diative cooling. 

This process produces a transient struc- 
ture, but this is not in conflict with esti- 
mates of only a few million years for the 
ages of the stars in Giant Molecular Clouds 
(GMCs) (Ballesteros-Paredes, et al. 1999; 
Hartman 2003). This is also not very dif- 
ferent from the time that a clump mov- 
ing at the observed velocities of ~ 10 km 
s _1 takes to travel a significant fraction of 
the size of a GMC (~ 100 pc). We shall 
therefore assume that GMCs are transient 
and that their formation is triggered by 
external disturbances. We make no at- 
tempt to follow the entire evolution of a 
GMC; instead we describe an idealised nu- 
merical calculation that illustrates how a 
cloud might make the transition from a 
state with moderate (3 to one with low (3. 

A significant fraction of interstellar 
space is filled by 10 6 K gas, which is 
composed of overlapping supernova rem- 
nants (Cox & Smith 1974; McKee & Os- 
triker 1977). This gas is continually be- 
ing reheated by shocks from supernova 
explosions. One estimate (McKee & Os- 
triker 1977) indicates that, on average, this 
causes the pressure at any given point to 
fluctuate significantly on timescales of the 
order of a million years or less. Such pres- 
sure fluctuations also appear in simulations 
of the structure of the interstellar medium 
(e.g. Gazol-Patino & Passot 1999; Korpi et 
al. 1999; de Avillez & Breitschwerdt 2004, 
2005; Mac-Low et al. 2005). de Avillez 
& Breitschwerdt (2005) and Mac-Low et 
al. (2005) do, indeed, find low (3 regions 
in their simulations. However, these are 
global simulations containing a number of 
ingredients, so that it is not entirely clear 
how these regions are produced. Our ide- 



alised calculation not only has the advan- 
tage that the mechanism is clear, but also 
makes it possible to use higher resolution 
than is possible in a global simulation. As 
we shall see later, even this higher reso- 
lution is not sufficient to follow the entire 
evolution accurately. 

2. The Model 

2.1. Heating and Cooling 

For the heating and cooling rates appro- 
priate for diffuse atomic gas, there are two 
thermally stable phases which can be in 
pressure equilibrium with each other for a 
range of thermal pressures. The low den- 
sity warm phase is at ~ 10 4 K, while the 
high density cold phase is at ~ 10 2 K. If 
molecules form, as they do in the clumps 
within GMCs, then the temperature of the 
cold stable phase decreases. The thermal 
pressure in the warm phase has a maxi- 
mum above which only the cold phase can 
be in thermal equilibrium. 

As a simple model of the thermal be- 
haviour of the material in the cloud we 
adopt the following form for the net heat- 
ing rate per unit volume 



H = p[0.015 - pA(T)] erg cm" 3 s" 1 



where, 






A(T) = 


3.564 10 16 T 2 - 21 


< T < 141 K 


A(T) = 


9.1 10 18 T 


141 < T < 313 K 


A(T) = 


1.14 io 20 T - 56 


313 < T < 6102 K 


A(T) = 


1.924T 3 - 67 


6102 < T < 10 5 K 


A(T) = 


1.362T- - 5 


T > 10 5 K 



(Sanchez-Salcedo et. al 2002). H is of the 
same form as that appropriate for the dif- 
fuse atomic gas (Wolfire et al. 1995). If 
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we ignore the dependence of the ionization 
fraction on temperature, then the pressure 
is 



V 



pkT 



m 



where m = 2 1CT 24 gm is the mean par- 
ticle mass. For this model, the maximum 
density and pressure of the warm phase in 
thermal equilibium are 0.5 cm -3 and 3051/c 
respectively. Here k is Boltzmann's con- 
stant. 

2.2. Initial Conditions 

It is well known that a transition from 
the warm to the cold phase can be trig- 
gered by a relatively small pressure in- 
crease caused either by flow convergence 
(Hennebelle & Perault 1999, 2000) or a 
weak shock (Koyama & Inutsuka 2000, 
2002, 2004). If is initially of order unity, 
then this transition can produce a small 
value of (5. We therefore consider a cloud 
which is initially in the warm phase and 
in pressure equilibrium with the surround- 
ing hotter material. The hot gas can be 
assumed to be adiabatic, but both heat- 
ing and radiative cooling are important 
within the cloud. We then suppose that 
the cloud is subjected to an increase in the 
pressure of the surrounding hot material 
from a value at which the warm phase ex- 
ists to one at which it does not. The pres- 
sure increase is assumed to take place on 
a timescale that is short compared to the 
time required for the cloud to adjust. 

In all the calculations, the cloud is ini- 
tially in thermal equilibrium in the warm 
phase with a number density of 0.45 cm -3 , 
corresponding to a temperature of 6277 K 
and a pressure of 2825/c. The external pres- 



sure is imposed by embedding the cloud 
in uniform hot material with density 0.01 
cm -3 with the appropriate pressure. This 
material is not in thermal equilibrium, but 
its cooling time is so long that it can be 
assumed to behave adiabatically. The ini- 
tial magnetic field is in the axial direction 
and is uniform everywhere with a pres- 
sure equal to the thermal pressure in the 
cloud i.e. P — 1 in the cloud. Since the 
cloud is initially spherical, the flow remains 
axisymmetric. Table 1 lists the external 
pressure and size of the cloud for the five 
cases considered. Changing the cloud ra- 
dius only affects the ratio of cloud radius 
to the cooling length behind the fast mode 
shock that propagates into the cloud. This 
ratio is ~ 0.1 for a cloud radius of 200 pc. 
Note that the external thermal pressures 
are well within the range found in global 
simulations. 

Model E explores the effect of an ex- 
ternal thermal pressure, p e , that decreases 
with time according to 



Pe 



p e */* c for p e > p 
Po otherwise. 



Here po = 10000A; is the initial external 
thermal pressure and t c = 2.6 10 7 yrs, 



External Gas Pressure 
Model A 15000A; 
Model B 10000A; 
Model C 10000A; 
Model D 7500A; 

(decreasing external pressure) 
Model E 10000A; 200 pc 



Cloud Radius 
200 pc 
200 pc 
50 pc 
200 pc 



Table 1: External pressure and cloud ra- 
dius. 



3 



which is approximately twice the time it 
takes the fast shock to return to the edge 
of the cloud (see section 2.4) . 

2.3. Computational Details 

All the calculations were done in ax- 
isymmetry using the unstructured hierar- 
chical adaptive grid code described in Falle 
& Giddings 1993. This uses a hierarchy of 
grids G° ■ ■ ■ G N such that the mesh spac- 
ing on grid G n is Axo/2 n . Grids G° and G 1 
cover the whole domain, but the finer grids 
only exist where they are needed. The so- 
lution at each position is calculated on all 
grids that exist there and the difference be- 
tween these solutions is used to control re- 
finement. In order to ensure Courant num- 
ber matching at the boundaries between 
coarse and fine grids, the timestep on grid 
G n is At /2 n where At is the timestep 
on G°. Such a grid structure is very effi- 
cient for flows that contain very thin dense 
regions since it confines the fine grids to 
where they are needed. The basic algo- 
rithm is an MHD version of the second 
order Godunov scheme described in Falle 
(1991). This uses a linear Riemann solver 
and the divergence cleaning algorithm de- 
scribed in Dedner et al. (2002). 

The computational domain is < r < 
2R C , < z < 2R C where R c is the cloud 
radius. The boundary conditions are sym- 
metry on the z = plane and the axis and 
free flow conditions elsewhere. The low 
density in the external gas ensures that the 
boundaries at r = 2R C and z = 2R C have 
a neglible effect on the solution. There 
were 6 grid levels with the finest being 
1280 x 1280. This gives a mesh spacing 
of 0.31 pc for models A, C, D, E and 0.08 
pc for model B. Even though this is signifi- 



cantly higher resolution than that attained 
by de Avillez & Breitschwerdt (2005) and 
Mac-Low et al. (2005) it is not sufficient to 
follow the complete evolution of the cloud. 
It is, however, enough to show the forma- 
tion of the fragments that we identify with 
the translucent clumps and to give some 
insight into their initial evolution. 

2.4. Results 

Figure 1 shows grey-scale maps of the 
number density, magnetic field lines and (3 
for model A at the times shown. A fast 
mode shock is driven into the cloud, which 
reflects off the symmetry axis and returns 
to the edge of the cloud at the last time 
shown in Figure 1 (13.2 10 6 yrs). At this 
point most of the cloud has (3 ~ 0.03 and 
is in the thermally unstable temperature 
range. Note that self-gravity can only have 
a minor effect upon the dynamics up to this 
point since even at the latest time shown in 
the Figure, the free-fall time is ~ 3 10 7 yrs, 
whereas the crossing time at the thermal 
sound speed is ~ 2 10 7 yrs and that at the 
fast mode speed is ~ 6 10 6 yrs. 

At this point the mean density of the 
cloud is ~ 30 cm -3 , which is typical of that 
in a GMC. Furthermore, the low value of (3 
provides the ideal conditions for the gener- 
ation of density inhomogeneities by magne- 
tohydro dynamic waves (Falle & Hartquist 
2002), particularly since most of the gas 
is in the thermally unstable temperature 
range. The end result must be the forma- 
tion of dense regions that are in the cold 
phase embedded in diffuse warm phase. In- 
deed, Figure 1 shows that such dense re- 
gions have already begun to form at the 
latest time shown. Unfortunately, we can- 
not follow the evolution of these dense re- 
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gions much beyond this because not only 
does the numerical resolution become in- 
adequate, but self-gravity begins to influ- 
ence the large scale dynamics of the cloud. 
However, it is clear that these regions will 
eventually come into approximate pressure 
equilibrium with the rest of the cloud, in 
which case their densities must be in the 
range 500 — 1000 cm' 3 . This is typical of 
the densities in the translucent clumps in 
GMCs (e.g. Williams, et al. 1995). 

There is also a slow mode shock which 
moves much slowly than the fast mode 
shock because the reduction in the mag- 
netic pressure across it means there is 
nothing to stop the gas being compressed 
as it cools, so that it comes into thermal 
equilibrium in the cold phase. This shock 
therefore produces a thin dense shell near 
the boundary of the cloud in which /3 ~ 1. 
Since the post-shock gas is in the ther- 
mally unstable regime, this process is un- 
stable and we would expect the thin shell 
to fragment. Our calculations show that 
this does indeed occur at later times. 

As one would expect, the increase in 
magnetic field inhibits collapse perpendic- 
ular to the field so that the cloud becomes 
flattened along the field lines. The slow 
shock must eventually eliminate the dif- 
ference between the thermal pressures in 
the surroundings and the cloud so that it 
ends up as a thin disc in which (5 is of 
the same order as that in the surroundings. 
However, this process is not complete until 
t ~ 2.4 10 7 yrs. 

The behaviour of the other simulations 
is similar, except that larger external pres- 
sures lead to faster evolution and smaller 
values of (3. In order to quantify this, in 
Figure 2 we show the fraction of the vol- 



ume of the cloud for which /3 < 0.1. It 
is clear that in all cases there is an ex- 
tended period during which (3 < 0.1 in 
much of the volume of the cloud. Figure 
3 shows that for models A and B, there is 
also a significant fraction of each cloud for 
which /3 < 0.05. The results for model E 
show that there is still a reasonable volume 
of low (3 material even when the external 
thermal pressure decreases with time, al- 
though the timescale for this decrease must 
not be much smaller than that in model E. 
For example, there is little low f3 material 
when t c is reduced to half this value. 

Although a significant fraction of the 
volume of each cloud is occupied by ma- 
terial with low f3, the corresponding mass 
fraction is somewhat lower. This is be- 
cause, once the fast shock has returned to 
the cloud boundary, the cloud consists of 
three regions: the region interior to the 
slow shock, in which f3 is low and the gas is 
in the thermally unstable phase; the thin 
shell behind the slow shock in which /3 ~ 1 
and the gas is in the cold phase; a thin an- 
nulus on the symmetry plane in which the 
conditions are similar to those in the thin 
shell. Although the thin shell and the an- 
nulus contain a significant fraction of the 
mass, they only occupy a small fraction of 
the volume. 

The annulus is not very significant since 
it is an artifact of the symmetry of the 
calculation and would not be present in 
a more realistic calculation without this 
symmetry. On the other hand, a thin shell 
must appear in any such calculation. This 
is interesting since it is likely to break up 
into dense fragments (Koyama & Inutsuka 
2000) and is therefore a possible site of 
star formation. In a more realistic simu- 
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lation which is asymmetric due to an ex- 
ternal flow, it would not be as regular as 
in our simulation and so one might ex- 
pect star formation to occur on one side of 
the cloud. There is a good deal of obser- 
vational evidence that suggests that high- 
mass stars tend to form at the surfaces of 
GMCs ((e. g. Israel 1987; Gatley et al. 
1979: Fich et al. 1982). In the case of 
the Galaxy, this evidence is not conclusive, 
but observations of other galaxies may pro- 
vide clearer results. For example, there 
is some tentative evidence from M33 that 
young high-mass stars are slightly off-set 
from the molecular material (Pattison & 
Hoare 2005). 

3. Conclusions 

The above simulations show the for- 
mation of highly magnetically dominated 
regions in thermally unstable material be- 
hind fast-mode shocks propagating through 
a cloud that is initially in a warm, ther- 
mally stable state. The external pressure 
variation required to generate such regions 
is well within the range expected in the 
mid-plane of the Galaxy. A slow-mode 
shock follows the fast-mode shock lead- 
ing to the production of dense regions in 
which the thermal and magnetic pressures 
are comparable. These dense regions may 
be the sites of high-mass star formation. 

We thank L. Blitz and M. Hoare for in- 
put and the referee for helpful comments 
on the original version. A. J. Lim ac- 
knowledges support from PPARC during 
the course of this work. 
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Figure 1: Linear grey scale plots of the 
number density and magnetic field lines 
(left) and f3 (right) at the times shown. For 
the number density the range is - 4 and 
for it is - 1. 
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Figure 2: Fraction of the volume of the 
cloud for which (3 < 0.1. 
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Figure 3: Fraction of the volume of the 
cloud for which (3 < 0.05. 
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